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A CLASS OF PERIODIC ORBITS OF AN INFINITESIMAL BODY 
SUBJECT TO THE ATTRACTION OF n FINITE BODIES* 

BY 

WILLIAM EAYMOND LONGLEY 
Introduction. 

The question of periodic solutions of a set of differential equations with appli- 
cation to the problem of three bodies has been treated by PoiNCARfi,t and the 
idea of analytic continuation which he developed will be employed in this paper 
to prove the existence of the solutions under consideration. 

Suppose that n finite bodies move in one plane and that their coordinates 
with reference to one of the bodies, M, as origin are known functions of the 
time. At any epoch t^ the bodies form a certain geometrical figure and have 
certain relative components of velocity. Suppose that at the epoch t^ -f- T^ the 
bodies form again the same geometrical figure and have the same relative com- 
ponents of velocity ; it is said that the motion is periodic with the period T^ . 

A particle T* is supposed to move in the plane of the motion of the finite 
bodies, subject to their newtonian attraction. In the differential equations of 
motion of the particle it is convenient to divide the terms into two classes : 
(1) those which are due to the attraction of M alone, and, (2) those which are 
due to the attraction of the other bodies. A parameter fi is introduced as a 
factor of all the terms of the second class, so that when /a = the differential 
equations represent the motion of P when subject to the attraction of only one 
body, M. In this case the solution is known and the conic in which the particle 
moves is called the undisturbed orbit. 

If there exists a periodic solution of the differential equations when fi ^ ,X 
the period T oi the solution must be a multiple of T^, T= qT^. Suppose the 
period is assigned in advance. When ^= a periodic solution is known. The 
initial conditions may be determined so that the period T„ of this solution is a 
submultiple of the assigned period, that is, p2\— qT^= T, where p and q are 

* Presented to the Society April 14, 1906. Received for publication June 15, 1906. 

t Bulletin Astr'onomlque, vol. 1 (1884), p. 65 ; Acta Mathematica, vol. 13 (1890), 
p. 5 ; Les methodes nouvelles de la mecanique celeste, vol. 1, chap. 3 (1892). 

X The analysis employed in this paper is valid for sufficiently small values of the parameters 
involved. All statements concerning the properties of the solutions for /i 4= refer only to values 
of i« so small that the solutions are known to be valid. 
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integers. There exists, therefore, a solution o£ the differential equations for 
fj, = 0, which has the assigned period. Does this solution persist, in the sense 
of analytic continuation, for values of fi different from zero but sufficiently 
small ? This question is here investigated by the method of Poincae6. Three 
cases present themselves for consideration : 

I. When the orbit of the particle is symmetrical with respect to a line. 

II. When there exists a uniform integral of the differential equations of 
motion. 

III. When no uniform integral exists. 

I. The first case may occur as a particular sub-case of either the second or 
the third ; but, if the differential equations admit a symmetrical solution, the 
proof of the existence of a symmetrical periodic orbit does not depend upon the 
consideration of the integral, and may be most conveniently treated as a distinct 
case. 

If a particle is subject to the attraction of two finite bodies, there are always 
initial conditions under which the orbit is symmetrical with respect to the line 
joining the finite bodies, and with respect to the time of crossing this line. If 
the body J/j revolves about Jf in a circle, the condition for a symmetrical orbit 
of the particle JP is that P shall cross the radius vector of M^ orthogonally at 
any time. The orbits of Hill,* DAEWiN,f and MouLTON J are of this type. 
If the orbit of J/j is an ellipse, the condition is that P shall cross the radius 
vector orthogonally when J/, is at an apse.§ The property of the symmetry of 
the orbit of J* with respect to a line rotating with an angular velocity which is 
a function of f, holds for some cases in which the particle is subject to the 
attraction of more than two finite bodies. When the motion of n finite bodies 
is such that they form a geometrical figure having always the same shape, in 
which equal masses are symmetrically situated with respect to a line of the 
figure, the particle may be started so that its orbit will be symmetrical.! An 
example of another type of the motion of the finite bodies, which is such that 
the orbit of the particle may be symmetrical, is given in the last section of this 
paper in the third numerical problem. 

II. The case where a uniform integral exists occurs when the finite bodies 
move in circles about Jf with the same angular velocity W. The bodies form 
a fixed geometrical configuration in a plane which rotates about the origin with 

♦American Journal of Mathematics, vol. 1 (1878), p. 245. 

tActa Mathematical vol. 21 (1897), p. 99. 

tTransactions of the American Mathematical Society, vol. 7 (1906), p. 537. 

§ This case has been treated in an unpublished paper by W. D. MacMillan presented to the 
Society, April 14, 1906. 

II Examples, are furnished by the equilateral triangular solution of Lagrange, when two of 
the bodies are equal, and by many particular solutions of the problem of n bodies. See W. R. 
LoNOLEY, Some particular solutions in the problem of n bodies, Bulletin of the Amprican 
Mathematical Society, vol. 13 (1907), p. 324. 
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the angular velocity iV. When the coordinates of the particle are referred to 
axes in the rotating plane, the differential equations of motion do not contain 
the time explicitly, and the period of the solution may be chosen arbitrarily ex- 
cept for limitations which are necessary to assure the convergence of the series 
employed. 

III. The third case occurs when the finite bodies do not move in circular 
orbits about M. The parameter /i is introduced into the differential equations 
in such a way that the proof of the existence of periodic orbits of the particle is 
made to depend upon the terms which are due to the attraction of only one of 
the perturbing bodies, M^ . The motion of M^ is subject to mild restrictions, 
namely, that the radius vector shall be given by an expression involving only 
cosines of multiples of the time, and that the expressions for the longitude shall 
involve only sines of multiples of the time. The only restriction upon the 
motion of the other bodies is that it shall be periodic. The treatment is 
sufficiently general to permit applications in the problems presented by the 
motions of the solar system. For example, suppose that P is a satellite of 
one of the planets M, and that M^ is the sun. The conditions upon the 
motion of J/j are fulfilled if we neglect the perturbations of the other planets 
upon M. If we neglect the inclinations of the orbits of the other planets, and 
suppose that their motion is periodic, it is possible by the methods given below 
to treat the periodic motion of the satellite in the plane of the planetary orbit 
when subject to the attraction of the sun and all the planets. 

It is shown that when the period of the motion of the particle is assigned, 
there exist two and only two real orbits of this analytic type having the pre- 
scribed period. In one of the orbits the motion is direct, and in the other it is 
retrograde. If symmetrical orbits exist, there are no unsymmetrical periodic 
orbits. It is shown that, in the analysis employed, the eccentricity of the undis- 
turbed orbit must be zero. A method is given for constructing the periodic 
solutions to any desired degree of accuracy by processes which involve only 
algebraic computation. 

For conve lience and clearness the details of the analysis are carried out for a 
special example of the motion of the finite bodies, namely, for three finite bodies 
which move according to the equilateral triangular solution of Lagrange. This 
example illustrates the three cases which may occur, and only slight changes in 
the details are necessary in order to apply the treatment to the more general 
examples previously mentioned. 

Some numerical examples are given at the end of this paper with drawings 
of the orbits which have been computed. 
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Existence of Periodie Orbits. 

The differential equations of motion. Let M, M^, and M^ denote the 
masses of three finite bodies moving according to the lagrangian equilateral 
triangular solution. With reference to Jif as origin and an axis having a fixed 
direction in space, let the coordinates of M^, M^, and the particle P be respec- 
tively {Hi, T^i), (-ffj T^), and (r,v). The coordinates of the bodies are 
expressed in terms of the time as follows : 

Ii^ = Ii^ = AS = a\ l-ecosm-^^{cos2m-l) + ■■■], 

(1) 

F, — TBTj = Fj — -BTj = F= JVi + 2e sin M + |e' sin 2M + ■■■, 

where 

M'A' = l^{M+ Ml + M^). 

The finite bodies are at apses of their respective orbits at < = ; JVis the mean 
angular motion, A the major semi-axis, e the eccentricity, ^ is a constant depend- 
ing upon the units employed ; and w^, nr^ are constants such that vr^ — Wi = irjZ. 
The differential equations of motion of P are 



d^r /dv\ 

'df~''\Jt) 
(2) 



2 ¥M an 



d^v dr dv 1 dQ, 

''W^ dt di'^r ~d^' 



where 



,SMi M^ Ml M^ , __.1 

L^ ^ ~ :^'-cos(« - FJ --^'/-cos (^> - FJ J , 



(3) ri=Vr'-\-R\-1rRiCa?.{v-Vi), 



r^ = Vr^ + Pl- 2rP^ cos (u - FJ . 
Let us define m and a by the relations 
(4) mv = jy, i;='a' = FJf, 

where v is the mean angular velocity of P in the undisturbed orbit. 

The period of the motion of the particle in the undisturbed orbit is 27rfv. 
The period of the motion of the finite bodies is lirfN. In order that the periods 
shall be commensurable, v is determined by the relation 2p'jr/v = ^qirjN, 
where p and q are integers. This restriction is not necessary, however, when 
e = 0. 

The next step in obtaining the final form of the differential equations is 
to refer the motion to a plane rotating with the angular velocity of the finite 
bodies. If e = , the coordinates of J/J and J£, become constants, and v may be 
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chosen arbitrarily. The motion is referred to an axis rotating with the angular 
velocity of the finite bodies by the substitution v = w-{- V. It is convenient 
to eliminate the factors depending upon the units employed by, the substitution 

r = ap, vt=r. 

Hence Ni becomes ttit. We may expand li as a power series in apj Alt as 
follows : 

^^\Tr) {3cos(«)-Tirj)4-5cos3(M)-T!rj|+ . . . 1 

+ H^i) {^°«K^-'^2) + 5cos3(w-w,)j- + ...J 

Let \ and \ be defined by the relations 

(5) M, = X,{M, + M,), M, = \{M, + M,). 

From equations (1) we have 

]c\M^ + M^) N^ 

A' ~ M • 

On setting 

-K, 



1+ ^ 



it follows that 



+ -^' (^) { 3 ''o^ (M,--nr,) + 5 cos 3(w --etjI + 



(6) 



+ -ii?{H-3cos2(w--Brj| 



■•■ 8 VZ&){^°°^^'"~'^2) + 5cos3(M)--or,)| + •• j, 
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1 an Km'prs\ 

I Ysin2(w-i!rJ 



(6) 



a^p dw R^ 

+ -Y-sin2(«)--nrJ 



The differential equations of motion become 

"5? ~ '' V c^T "^ dT ) '^ p^~~a F{ap) ' 



/^ ^F\ dpfdw dV\_ 1 aft 

'^ V'^ + ^7V '^ rf^ \^ '^^ )~ aV awr • 

Case /. Symmetrical solutions. Suppose the masses iH/^ and J/^ are equal 
and that the w-axis passes through the center of gravity of the system. These 
conditions are expressed by the relations, 

(8) A., = \, tD-j = - - , ^2=6- 

It follows that the first of equations (7) involves only cosines of multiples of w, 
and the second involves only sines of multiples of w . Since d V/dr and li are 
even functions of t, it follows that the equations are unchanged if p is left 
unchanged while w and t are changed into — w and — t respectively.* When 
subject to the restrictions (8), let the differential equations be denoted by (Ti).t 
Suppose that 

(9) p = yjr^{T), w = ^|r^{T), 

is a solution of (Tj) such that 

(10) /- =0, w{pTr)-p7r=0. 

The differential equations are periodic in w and t, and JJ/j and J/^ are at apses 
at T = pir . It follows from the form of the differential equations that ^jr^ is an 
even function and yjr^ is an odd function of t . Hence if the particle crosses the 
to-flxis when the finite bodies are at apses of their respective orbits, the orbit of 
R is symmetrical with respect to the axis and the time of crossing. If the 
orbit is symmetrical with respect to two epochs, for example, t^^O and T,=j97r, 
it is periodic with the period 2297r. If the solution (9) is symmetrical with 

* In the more general cases in which symmetrical solutions occur it may happen that the first 
equation contains also sines of multiples of w multiplied by odd functions of t, and that the 
second equation contains also cosines of multiples of w multiplied by odd functions of t . 

t Throughout this paper the subscripts I, II, III refer to the wise I, II, III, respectively. 
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respect to the epoch t = , the conditioivs (10) are necessary and sufficient for it 
to be periodic with the period 2pTr . 

We will now investigate the existence of symmetrical periodic orbits which 
have the period 2/>7r. The terms of equations (Tj) are periodic in t and ana- 
lytic in m ; but, when expanded in this parameter, the coefficients of the various 
powers of m contain non-periodic terms. It is convenient in this discussion of 
periodic solutions to avoid such terms in the differential equations, and this is 
accomplished by a generalization of the parameter m.* Consider the set of 
differential equations, 



(11) 
where 



d?p (dw duy 1 

fd'w d^U\ ^dpfdw dU\ 
P\d?-^''l[^)^^dr\d^^''-d^) = ''^^^ 



r/= — V= — j mr -f- 2e sin nn + i^ sin 2mr -f 
m 7n '- * 



and /x^y*! and p^g^ are obtained from \ j a ■ dCl j d {ap) &i\A l/a/'p dil/dw respec- 
tively by writing p,^ for the factor m\ and by making the substitution a/^l = rjp,, 
where t; is a numerical constant. 

Equations (11) do not represent any physical problem except for a sjiecial 
value of p, namely, p = m, and are to be studied from a mathematical point 
of view. They may be solved for p — 1 and w as power series in p. 
For a preassigned arbitrary interval for t the series are convergent for all 
values of p numerically less than /i,| and satisfy the differential equations 
identically in p. Therefore, if in the solution p is given a value numerically 
equal to m(^m<:Cp), the solution satisfies the differential equations (7,) and 
represents the physical problem under consideration. It follows from the form 
in which m enters the differential equations that p is independent of ?n . 

When p= the equations are of the form occurring in the problem of two 
bodies, and a symmetrical solution having the required period is known, namely, 

p = 1 — e cos B, 

/ cos ^ - e \ . / v/r^e^ sin ^\ t 

I ., — ^ I = arc sin I ^ v^- I, 

\ 1 — e cos ^/ \ 1 — e cos i / 

where £J is defined by the relation 



(12) / cos ^ - e \ . / V^l - e^ sin E 

w = arc cos 



T = ^ — e sin ^, 



*F. R. MouLTON, On certain rigorous methods of treating problems in celestial mechanics. 
The Decennial Publications of The University of Chicago, vol. 3(1902), p. 126. 
tPoiNCAEE, Lea methodes nouvelles de la mecanique celeste, vol. 1 (1892), p. 58. 
t The positive sign is to be taken with the radical. 
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and e is an arbitrary constant less than unity. The initial conditions for t = 
are 

(13) p=l-e, ^^^ = 0, ,0 = 0, ^=^^_^. 

C!onsider the solution for values of /t different from zero but sufficiently 
small, and let the initial conditions for t = be 

p=l-e + ^ = (l + a)[l-(e+e)],* J = 0. ■"' = 0. 

(14) 

dw _ V\ - e' _ v/l-(g + e)' dD^ 

dr~{l-ey'^'^~ [l_(e + c)]^~'* dr 

The solution is symmetrical with respect to the epoch t = , and may be 
expressed as power series in a , e , and /x , which are convergent for an interval 
in T including the interval to pir , if the parameters are sufficiently small. If 
a and e can be determined in terms of fj., vanishing with fi, so that the condi- 
tions (10) are satisfied, then the solution will be periodic with the period ipir. 
All terms of the solution which are independent of fi^ may be obtained from the 
two body problem by making the substitution to = u ~ nU. 

These terms are given in finite form by the expressions 

p==(l + a) -(e + e)cos J?] , 



%i = arc cos 



. co8^-(e + e) \ . ( \/l-{e + ef9,\nE \ 

\l-{e^e)cosE)~'^^^^^'^\ 1 - (e + e) cos jS" )' 

where E is defined by the relation 

Returning to the variable m), writing the terms in a and e as power series by 
Taylor's expansion, and applying the conditions (10), we obtain the equations 

rt o V 1 — e'' 

= - f p,r __-___^.^a _^^^ + . . .. 

It follows from the known properties of the series that there are no terms in e 
alone, and there ai-e no terms involving fi to the first degree except the term 

*The introduction of the very eoovenient parameters a and «, instead of the additive incre- 
ments /J and y , is dne to Professor P. R. Moulton, these Transactions, loo. cit., from whom 
I have received many suggestions and much valuable criticism. 
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— pirii. Hence the second of equations (15) may be solved for a as a power 
series in e and fa in which /u is contained as a factor : 



(16) a=M|- 



(1 - ef 



When this value of a is substituted in the first equation, a factor fi may be 
divided out, leaving 

This equation may be solved, for e as a power series in ft, which vanishes with /*, 
if and only if e=0, and the result may then be used to eliminate e from equation 
(16). Since only those solutions are under consideration which are the ana- 
lytic continuations with respect to m of those for ^ = , the condition e = 
must be imposed. The condition e = means that the undisturbed orbit must 
be circular. 

The initial conditions (14) have been determined uniquely as power series 
in fi, convergent for sufficiently small values of /it, so that the conditions of 
periodicity are satisfied. When the values of a and e thus found are sub- 
stituted in the solution, we have p — 1 and w represented as power series in n 
alone, and the solution is periodic with the period ipir. 

Suppose /u = TO and consider the physical interpretation of the solution. The 
period of the motion of the finite bodies is 2'jr/m = 2pTr/q , where p and q are 
integers which may be chosen arbitrarily except for the condition that q/p must 
be sufficiently small. The period of the solution is 2pTr. Hence during this 
period the finite bodies make q revolutions, the angle w is increased by 2p7r , and 
therefore the particle makes q + p revolutions with reference to fixed axes. If p 
and q have the same signs, that is, if m is positive, the motion of the particle is 
direct. If »i is negative, the motion is retrograde. For a given value of to 
there exists one and only one real symmetrical periodic solution of the differ- 
ential equations. Hence for a given period there exist two and only two real 
symmetrical orbits of the particle with the required period. In one the motion 
is direct, and in the other it is retrograde. All distinct orbits are obtained 
when p and q are prime to each other. 

The initial conditions were limited so that the solutions should be symmet- 
rical. When this condition is not imposed, and whether or not the differential 
equations admit symmetrical solutions, it will be shown in the next two sections 
that, for a preassigned period, there exist two and only two real periodic orbits, 
in one of which the motion is direct and in the other retrograde. Hence, if 
symmetrical periodic orbits exist, there are no unsymmetrical periodic orbits of 
this type. 

Trans. Am. Math. Soc. 18 
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The existence proof was made for a particular example, but did not depend 
upon the right hand members of the differential equations and consequently 
holds for any case in which the orbits of the perturbing bodies are such that the 
differential equations admit symmetrical solutions. The only variations which 
it is necessary to make in the details for the other cases, are in the particular 
manner in which the parameter fi is introduced, and in the choice of the rota- 
ting axis. It may happen in some cases that there exist several possible axes of 
symmetry. 

Case II. When a uniform integral exists. Suppose that M^ and M^ are not 
in general equal and that they move in circles about Jff with the ?^axis passing 
through M[ . These «onditions are expressed by the relations 

(17) 6=0, •Br, = 0, •=^, = J. 

When subject to the conditions (17) let the differential equations of motion be 
denoted by ( Tn ) . 

In this case it is convenient to introduce the parameter /i by the relations 

a 

(18) m = fi, X^ = \/j,, -^ = rjn, 

for all powers of a/A higher than the Jirst. Here it is not necessary to gen- 
eralize the parameter wi since H and dVjdr are constants. By relating \ 
and n, the existence proof is made to depend only upon those terms of the dif- 
ferential equations which involve \. The generalization of the parameter ajA 
is merely for convenience in having finite expressions in the equations which 
determine the coefficients at the various steps in the solution. 
The differential equations become 



(19) 



where 



<Pp (dw X 1 
d^w ^dp (dw \ . 



/ii = ^p \| 1(1 + 3cos2w) + |.^/3(3cosw + 5cos3w) 

-\- ■■ ■ fi func (p, fi, coajw, sin jw) I , 

gji = — ICp X^\ I sin 2w + I ^ p (sin w + 5 sin 3io ) \ 

+ ■ ■ ■ fi func (p, fi, sin Jw, coajw) I . 
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Equations (19) are periodic in w with period 2ir and do not involve t 
explicitly. Suppose that 

p = -./r,(T), w = ir^{T), 

is a solution. The necessary and sufficient conditions that the solution shall be 
periodic with the period 2pir are 

,^,(2p.) = ^.(.0), f l^,^ = f [, 

(20) ' " 

When /t = a periodic solution is known,* namely, 

P = 1, M = T. 

The initial conditions for t = are 

1 '^P f\ n ^^ 1 

p = l, „ = 0, «, = 0, ^^^=1. 

For fjL different from zero, but sufficiently small, consider the solution of 
equations (19) subject to the initial conditions 

p = l + )3j = (l + a)(l-ecos^), 
dp e sin ^ 

<ir~ ^~ v'i +"a(i _ e cos 6) ' 

(21) ^ r cos0-e 1 [- cos(6-S)~e 1 
w = p, = arc cos :7 ^ — arc cos ^ ^ y^ rx , 

' LI — e cos » J LI — e cos (p — <^) J 

dr ^ * (l + a)i(l-ecos6l)^ 

The parameters, /3, , /S^, /S^, /8^, which are additive increments to the initial 
conditions of undisturbed motion, are inconvenient in the discussion which 
follows. In terms of the parameters a, e, 6 , 4>, the properties of the solution 
are well known, and the conditions of periodicity can be easily discussed. The 
geometrical meaning of the angles d and <f> is shown in Fig. 1. 

We may now write, in finite form, those terms of the solution which do not con- 
tain fj-, namely, 

r = (l +a)(l — ecos^), 

cos JE — e 

(22) cos{v, + W,-W,) = ^^^^^, 



^ ' ° -^ 1 — e cos M 



* As in the case of symmetrical orbits, no greater generality is obtained by assuming that the 
undisturbed orbit has an eccentricity diflerent from zero. Therefore we will assume here that 
e = 0. 
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where E is defined bv the relation 



(l + «)' 



+ ^ — e sin ^ = ^ — e sin JE". 



The only term of the solution which involves ix to the first degree is —rfi, 
which occurs in the expression for w . To find the terms in /u.^ and /i*(^ we may 
write 

/0=1 +p^fl^ + pfl'<f)+ ■■■, 

w = T — Tfi + (t> + Wj^fi^ + w/JL^tf H • 

Fig. 1. 




Po is the position of the particle at i" = 0. 
Wf, is the longitude of the particle at t = 0. 
Wi — Wd is the longitude of perihelion. 

On substituting these expressions in the differential equations, there results for 
the determination of p^ and w^ the following set of equations, 

-^^- - 2 -^' - 3/3, = -2"^ (1 + 3 cos 2t) + — g-i ^ (3 cos T + 5 cos 3t), 






3^\, . „ ZK\ a , . t • o X 

— o sin 2t Q— -J (sin T + 5 sm 6t). 
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It follows that 

+ if 2 '^ ^'^ '^ ~ cos 2t — If -J cos 3t , 

+ J/ ^ T cos T + -Y- sin 2t + Y/- ;i sin 3t . 
By a similar computation it is shown that 
p = K\^-(^i + 2J^6- ^ ^ sin T + if ^ T cos T + 2 sin 2t + if ^ sin 3t1 , 

— if- -7 T sin T + -Y- cos 2t + || -3 cos 3t . 

The terms independent of n are obtained from equations (22) by Taylor's 
expansion, and the solution becomes 

/o = 1 + a — e cos t — ae( cos t — |t sin t) + e^ sin t 

+ ae^( sin t — I-t cos t ) ■}- p^ti' + Pt^^<t> + • • • 1 
(23) 

w ^ T — |tcx + 2e sin t + (|) — r/i — Saer cos t — e^ ( 1 — 2 cos t) 

+ e<f) + ae6(5 cos t + 3t sin t) + w^/i^ + w/i'(j) + • • • . 

Applying the conditions (20) that the solution shall be periodic, we have 

= — Zpiraed + ^-f -3- pirp? <|) + • • • , (a) 

= ?,inrae + J/- ^ ^97^== + • • • , (6) 

= — ipira — 2^7r/i — Gpirae + • • • , (c) 

= Qpiraed — -^5. —pTj-fi^ ^ ^ . . ., (<;) 

The conditions (24) involve the four quantities a, e, 6, <f>, and, if indepen- 
dent, would determine them in terms of fi. But the differential equations (19) 
do not involve t explicitly and hence admit the integral of Jacobi. This fur- 
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nislies a relation of the type 

F{a, e, d, <}>, /^) = constant, 

and equations (24) are not independent.* It follows that if (a), (b), and (c) are 
solved for the three quantities a, e and in terms of /x and <^, and the results 
substituted in (d), the equation is satisfied identically in ^ . In this problem 
the dynamical interpretation is simple. Since the finite bodies move in circles 
the origin of time is arbitrary. The most convenient choice is t = when 
M) = , which is equivalent to choosing <^ = . 

Consider the solution of equations (a), (h), and (c) for a, e, and 0. The equa- 
tions have the following properties : 

i. There are no terms independent of a and /m. This follows from the fact 
that, in the two-body problem, the period does not depend upon e and . 

ii. There are no terms involving /i to the first degree except the one term 
— 2^7r/x which occurs in (c). 

iii. There are no terms in independent of e, since does not enter the 
initial conditions independently of e. 

It follows from these properties and the particular form of the first terms of 
the equations that a , e , and may be determined uniquely as power series in fi 
by the following steps: 

(1) From (c) we obtain 

« = /*[— f H |-func(M, e, 0)] . 

(2) This value of a when substituted in (6) permits a factor /* to be divided 
out. Then we may solve the result for e as a power series in /* and , in which 
/i is contained as a factor, 

(3) When the values of a and c are substituted in (a) a factor /jl^ may be 
divided out and obtained as a power series in /t alone, vanishing with /t . 

(4) By substitution of the value of thus found in the expressions for e and 
a, we obtain finally 

The preceding operations are known to be convergent for all values of a , e , ^ , 
and /i which are sufficiently small. Hence for a given value of /jl sufficiently 
small it is possible to determine the initial conditions (21) as power series in /jl 
such that the solution of the differential equations (19) shall be periodic in t 
with the period ipv. 

Since iMv = iV, M is determined by choosing the period of the solution in t 

*Ct. PoiNCARE, Les mithodes nouvelles de la micaniqtce celeste, vol. 1 (1892), p. 87. 
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which is ^pnrjv. All distinct orbits are obtained by taking p = \, that is, 
every orbit is reentrant after one revolution in the rotating plane. For an 
assigned period there is one and only one direct orbit, and one and only one 
retrograde orbit. 

The analysis and results apply to the case when any number of finite bodies 
revolve about M in circles with the same angular velocity. 

Case III. When no uniform integral exists. Suppose that M^ and J^ are 
not, in general, equal ; also, that the eccentricity of the orbits of the finite bodies 
is not zero, and that the w-axis passes through J/^. These conditions are expressed 
by the relations 



•JT 



w„ ^ 



(25) e+0, ■Br, = 0, ^,-^. 

When subject to the conditions (25), let the differential equations of motion 
be denoted by (Tm). In this case it is convenient to consider the equations 

(26) ^ ^ 






where 17= fi/m- F" and M^/in and fi^gm are obtained from l/a-5fl/5(ap) 
and 1/a^p ■ dil/dw respectively, by writing jjl^ for the factor m^ and making the 
substitution 

\ = X/i. 

Equations (20) express the conditions that a solution shall be periodic with 
the period 2pTr. Consider the solution of equations (26) subject to the initial 
conditions (21), except that the fourth condition contains the term —fi-dU/dr 
(t= 0) instead oi — fi. The terms independent of fi are given by equations 
(22). The only term involving fi to the first degree is — n/m ■ U, which occurs 
in the expression for w . 

For the determination of the terms in fi^ we have the equations 

|^_2'^^-3/,,= ^'rj(H-3cos2T)+f~(3cosT+5cos3T)+...], 
(27) •- -^ 



^'■'^2 , o'^Pi ^^ 



' + 2-P = 



i fsin2T + |^(sinT+6sin3T) H . 



dr' ^ dr ~ B' 

It \B & power series in e whose coefficients involve only cosines of multiples of 
im . It follows that the right hand member of the first and second equations 
respectively have the forms 

Ao+ Z)^?,aC0s(i7 + Am)T, Y^ Dgj^%m {g -Jf h7n)r , 

g, A g, h 
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where g and h are integers, g taking all positive values and zero, and h tak- 
ing positive and zero values such that hm < 1 . Since t enters the differential 
equations explicitly, we have m = qfp , by equation (4). A^^ ^ and Dg ^ are con- 
stants depending upon e and a/ A; they are infinite series which, for small 
values of the arguments, converge rapidly. 

The second equation of the set (27) may be integrated. The result is * 

(28) ^^-^2p, = c.-i:^cos(, + A.)r, 

where c^ is the constant of integration. 

On eliminating dwjdr between the first of equations (27) and equation (28), 
there results an equation of the form 

-^ -f />j = J„, + ^1, cos T H \- A,^ ^ COB {g-{-hm)T + ■••', 

whence f 

/3j = Cj cos T + Cj sin T+ J„ o+ji[i,oT sin tH- \- j37^t|^2Cos(5'+^w)t- • •, 

where Cj and Cj are constants of integration. 

Substituting the value of p^ in equation (28) we get 

-^= A,o — 2c3 sin T — J^.oT sin T + 1- ^^,4 cos {g -\- hm)T + • • • ; 

whence w^ is obtained by quadrature. The constants of integration are deter- 
mined by the initial conditions, 

dp dw„ - 

It is important to observe that all the terms of the solution of equations (27) 
have the period 2/)7r, except the terms t and tcost which occur in w^, and 
T sin T which occurs in p.^. Consequently when the conditions of periodicity 
(20) are applied, equations (a) and {d) will contain no term in p? alone ; (6) and 
(c) wiU each contain a term in fi^ alone. 

For the determination of the terms in p,^^ in the solution we have equations 

of the form 

d^p „ dw „_ v-« ^ • / 7 V 
d^'^Tr -3p=E<5^i7,*sm(gr + Am)T, 

(29) 

-^ + 1-£^H,^, + 23/f^^ cos {g + hn)r. 



* By the definition ot g, h, and m the denominator g + '»»> cannot vanish. 

t By the definition ot g, h, and m the denominator 1 — (g -\-hmy cannot vanish. 
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The solution of the set (29) is of the form 

— (^ 

P = 2^o,oT -JG'i.o'r cos T + . • . + ]^_ . y^^x2 sin {g + hm)T +■ ■ ., 

— ff 

w = - l^coT^* + (?i,„T sin T + g^'hm «>s (S' + '^»i)t + • • •• 

The conditions that the solution of equations (26) under consideration shall 
be periodic with the period 2p7r are 

= - Sp-n-aee + (4ffo,„P'^ - ^i,oPt)m'<^ + • • •, (a) 

= — Spirae + 1^ ^pirix^ + • • • , (6) 

(30) 

= — 3p7ra + 2p'7r/t + • • • , (c) 

= Qp-iraed - {6H,^ .pir -2Gi^„pTr)fi^<}) + ■ ■■. (d) 

Equations (30) have the following properties : 

i. There are no terms independent of a and /t . 

ii. The first power of /it does not occur in (a), (6), and (d). 

iii. There are no terms involving <f> which are independent of /jl^ . 

iv. There are no terms involving which are independent of e. 

V. Equations (a), (6), and (d) contain no term in a independent of e. 

vi. In equations (a) and (d) there are no terms of degree less than the third. 

The knowledge of the preceding properties leads to the conclusion that a, e, 0, 
and <f> may be determined uniquely as power series in /t, vanishing with /j- and 
satisfying equations (30). This conclusion is reached by the following steps 
which are valid for sufficiently small values oi a, e, 0, <f>, and fi. 

(1) Let equation (c) be solved for a as a power series in e, d, <f>, and fi 
(which contains /* as a factor by i), giving 

(31) « = |M + /t^.(e,^,<^,M). 

Equation (31) may be used to eliminate a from equations (a), (6), and (d). 
Let the resulting equations be denoted by (a), {i'), and {d') respectively. 

(2) Then (6') contains a factor fi which may be divided out, leaving 

= — 2p-7re + AifipTTiM + • • • . 

This equation contains no term independent of e and fi (by i, ii, and v), and 
contains a term of the first degree in e . Hence it may be solved for e as a 
power series in ^, <^, and fi, which vanishes with /t, 

(32) e = ^l,^of^ + ,i%(0,<t>,fi). 
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Equation (32) may be used to eliminate e from equations (a') and ((?'). Let 
the resulting equations be denoted by (a") and (<?") respectively. 

(3) Each of equations (a") and {d") contains fj,^ as a factor (by i, ii, iii, 
iv, V, vi), which may be divided out, leaving no term independent of , <f>, 
and fj. : 

= - pttA^ „e + {iff^^p-TT -Gi^op7r)4> + ..., (a") 

(33) 

The jacobian with respect to and <f> for = (f> ^ im = is— 22>^7r^-4, „^ o, 
which does not vanish identically in e and a/ A . Hence equations (33) may be 
solved uniquely for and ^ as power series in /t, and the result used to elimi- 
nate and (/> from equations (31) and (32). This yields the final solution of 
of equations (30), 

(34) « = m^,(m), e = /.qj,(M), = M-%M, «^ = M^,(M). 

When these values are substituted in the solution of equations (26), the ex- 
pressions for p and w become power series in /t alone, and the solution is periodic 
with the period 2p'7r, since the conditions (30) are satisfied. The physical inter- 
pretation is the same as in the preceding cases. 

For the purpose of constructing the solutions in applications, it is convenient 
to have the existence proof in the form given. It leaves open, however, one 
question which should be answered. That is, are there any values of e different 
from zero for which the jacobian of equations (33) vanishes ? The jacobian is 
a power series in e and a/A, vanishing with e, and it is not easy to discuss 
other values of e for which the series might vanish. Such a discussion is made 
unnecessary by a slight variation in the details of the existence proof. The 
proof depends only upon the properties of the solution of the problem of two 
bodies and the non-periodic terms which enter in p^, w^, p, and w, and may be 
carried tJirough without considering in equations (27) and (29) the infinite series 
in e which represents IfH and the infinite series in a/A . This is accomplished 
by picking out from /,„ in (26) a single term of the type 

/ a \''-' 
€''1-^1 cos pmr cos (q — l)w 

and from ffm & single term of the type 

1 
cos ^))«T sin {q — 1 ) w - 



ar 



These terms are left unchanged, while in all other terms e is replaced by <r/i and 
a/ A by vfi. The set of differential equations thus obtained may be treated by 
steps similar to those employed in the consideration of equations (26), and a set 
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of equations corresponding to (33) will be obtained. The jacobian of this set is 
e^'^ajA)-'''-, except for a constant factor, and can vanish only if e = . 

Construction of the Solutions. 

It has been demonstrated for each of the three sets of differential equations 
(11), (19), and (26) that, for a given value of /m sufficiently small, it is possible 
to determine the initial conditions uniquely as power series in /a such that the 
solution as power series in fi shall be periodic in t with the period %p7r . A 
method will now be given by which the solution to any desired number of terms 
may be conveniently constructed. It is not necessary to determine the initial 
conditions explicitly in advance and the computation involves only algebraic 
processes. It will be convenient for brevity in the treatment to consider first 

Case III. When no uniform integral exists. Consider equations (26), and 
let the periodic solution, which is known to exist, be written 

^ggv p-l = p,fi+ p^fi^ + p^fi' + ■■■+ pX + ..., 

W — T = W^fl + W^p? + W^fJ? + ■■■ + W^p.* + ■ •■. 

Since the solution exists and is periodic for all values of ft. sufficiently small, 
each coefficient is periodic with the required period. 

Let the solution (35) be substituted into the differential equations (26) and 
the result arranged as power series in fi . The terms of the left-hand members 
have the following forms. (The accents indicate derivatives with respect to t.) 

2 . 5 . 

U =T + U =T A — e sin mr + -. — e^ sin 2mT + • • • , 
TO 4m 



d'p 



^,=1 + ^', ^^-=u'\ 



^^2 = p'l'P' + p2>' + ^y + • • • + p^v + • • •, 

+ Ip^ + 2K + ^p,{w[ + 1 + ^') + («.; + 1+ D-'y-\ M^ 

(36) + |>3 + 2w',+ 2(7o; +l + U'yw',+2py„+2pJiw[ + l+ €") 

+ p.(M).+l+ Uy]p>'+- ■ ■ + [j>,+2w', + 2{w[+l+ ^>u 
+ 2p._^{w[ + 1 -t t/')+2py^_^+ ...]/i<+..., 
^ = 1 - 2p^^ -{2p,- ^p\)p? -{2p,- Gp,p, + ipl),.' 

+ ...+ {2p,-6p._^p^+...)p.'+.... 
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+ [«»;' + <p, + p2«' + ^")i P''+ ■■■ 

+ [<' + w^liPi + • • • + P,-,{w[' + C")-] /*• + •••, 
dp / dw dU\ , ^ , ,,, ^ fr,^-, , 

+ [/>;+p;(«'; + i + c^') + p>;] /*' + ••• 

+ ip'i + p'i-i{w[ + ! + £?") + ••• + /)>;_,] /*• + •■•• 

The equations hold identically in /i. Hence, equating the coefficients of the 
first power of /*, we have for the determination of p^ and w, the following set 
of equations : 

It follows that 

p, = 2 ( 1 + cV>) + d'> cos T + <'' sin t, 
(38) 

w, = c^' - ( 4 + 3c<,')) T - t?'- 20(2" sin t + 2c<3') cos t, 

where c'/\ c'^'^ c^^ and c'^'^ are constants of integration. Since p^ and w^ are 
periodic, the coefficient of t in Wj must vanish. This condition determines the 
constant c'j'\ namely d^^ = — ^. The other constants of integration are so far 
arbitrary. 

On equating the coefficients of the second power of p., the following set of 
equations is obtained : 

^P^_2^^-%p, = {w[ + l + U'f + 2p^{w[ + \+U')-^p\+U,o. 
(39) 

5^ + 2^^^ — p^(«,;' + ^")_2/,;(«,; + i + £?") + ^„,o, 

where yiii and j'm.o ^^re obtained from yju and g^^ respectively by writing 
/ir=0,w=T,/3=l. The right-hand members are known functions of t and 
the equations have the form 

%- - 2 J-^ - 3p, = Af\ + (-VH" + ^i%) cos T + -VH'^ sin r 



V ' rl'^ ■"< J/3, 



+ ■■• + .4(,%cos(5' + Ato)t-|- 

^ + 2 I~ = (-V^^" + A%) sin T - Y-c'a" cos r 

+ ...+Z><,%sin(g' + A»i)T + 
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It is to be noted that c^^ does not occur in equations (40). Integrating the 
second equation we have 

^^ + 2p, = cf - ( V-c(" + A!U cos T - l^^p sin T 

(41) "^^ ^., 

On eliminating dwjdr from the first of equations (40) by means of equation 
(41), there results 

^^ + P2= ^«!o + 24^' + [- 2cy> + Jl<fo - 2i)i%] cos T - 2cy) sin r 

(42) ^ 22/2) \ 

In order that the solution of equation (42) shall contain no non-periodic term, 
the coefficients of cos t and sin t must vanish, whence c'^'^ and c^' are determined 
by the conditions 

With these values of c'j'' and c^j", the solution becomes 

p^ = ^(o?Q + 2cf > + c^/) cos T + (*=> sin T + . • . + 4% cos (gr + Am) T + • • • , 

where 

1 



dP,. = 



X^'-" g + hm)' 



Substituting this value of p^ in equation (41) and integrating, we obtain for w^ 
a solution of the form 

M,^= 61^- (2^1,% + ^<^?)r - 2c'/' sin t + 2(4^' cos t + • • • + ^^\ sin (5r+ Am)T+ • ■ • , 

where 

«'.-^(^.+^<'.)- 

Since Wj is periodic, c'j^' is determined by the condition 

2^<,% + 3cf = 0. 

Of the eight constants of integration, which have been intioduced in the first 
two steps, four (c^/>, c'^'?, 4'', cf ) have been determined uniquely and four 
(c^^'-, c^/', c<3^>, c'/') are still arbitrary. 

By equating the coefficients of the third power of /i the following set of 
equations is obtained : 

^ - 2 ^ - 3p3 = 2(«,; + 1 + ^' + pO«'; + ^pA< + 1 + ^') 

(43) + p,{w[ + 1 + U'f-6p^p^_ + ipl +/„,,!, 

'^ + 2 ^ = - />,<-(«,;' + c^> - 2p;(.«; +i + u')- 2p>; + 5r,„, , , 
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where ^/m i and gm^ i denote the coefficients of fi in /"m and g^ respectively. 
Equations (43) do not contain c^p. The right-hand members contain (f^^ 
linearly. The equations have the form 

+ 1- ^fn COS iff + hm)T + Bf^,, sin (g + hm)T + • • • , 

-5-* + 2% = C^\ + (-V-cf + I>/^>o) sin T + ( - -V-cf + e?,'o) cos t 

-f- • • • D^^,. sin (jc + hm)T + C^/_',. cos (p- + Am)T + • • -. 

Here C'^% involves c'/' lineai-ly. To avoid non-periodic terms in the solution, C-po 
must vanish, and this condition determines c^^ uniquely, if the coefficient of c'^'^ 
does not vanish. The coefficient is a power series in e and a/ A which vanishes 
with € but not identically in e. The question concerning special values of 
€ for which it might vanish is the same as the question concerning the vanish- 
ing of the jacobian of equations (33) and iis answered in the same way. 

The treatment of equations (44) proceeds by steps similar to those employed 
in the solution of equations (40). Four new constants of integration are intro- 
duced, namely, c',^\ c^", di\ c'-^^ ; while c'^ ', d^\ and (^^^ are uniquely determined 
by the conditions : 

2cf = ^(■'.>o - 2J9'^)„ , 2cf = 5(,% -I- 2 C^,% , 3cf = - 2^«>„ . 

So far c^p, dp, cf^\ and c'^^^ are arbitrary. 

It may be established by complete induction that the above process can be 
carried as far as is desired. Suppose p,, w, ; Pj' ^2 ' ' ' Pi-i' '"'i-i ^^^® heea 
determined by this process. The expressions have the following form 

Pi = <'« + f- <\eos(gr-f ATO)T-]-/3(;_',.sin (gr-f Aj?i)T-f • • •, 

0=1, 2,- • -,1-3). 
Wi = 7'o?o + 1- oy,\ sm (g+hm)T+r/l>^ cos {g+hm)T+ • • •, 

Pi-2 = <■;;" + • • • + <'l'^ cos(g + hm)T + ^-^-iP sin (gr -f Am)T -I- • . ., 

w,_, = cy-« -I- • • • -t- r-^> sin {g + hm) r + i^;^ cos {g + hm)T +.-■. 

Pi-i = <i;" + c<;-'> cos T 4- 4'->> sin T -H • • • , 

w._j = c^;-'^ — 2cy-'> sin T -I- 2c*;-" cos t -( . 

The equations for the determination of p^ and w^ are 

2^ - ^ W - ^''*=^«?o+ (¥«';-"+ ^^?o ) cos T + ( -V-c(*-')+^,')„ ) sin T 



(45) 



+ ■■■ + Af;^„cos{g + hm)r -f R^^^^sm{g + hmjr + ■ ■ 

^i? + 4t = ^»- » + (-V-^^" + ■^'- « ) si'^ -^ + ( - -'M-'' + ^i?o ) cos r 
+ h 2)f,*>j sin (g + hm)T + C^^\ cos {g + hm) sin t -f- • • 



1907] SUBJECT TO THE ATTRACTION OP W FINITE BODIES 181 

Equations (45) do not contain c^^*~'\ The coefficient of c'^'~^' in Ci*jo is the same 
as the coefficient of c^' in C^^^o . Equations (46) are solved by the steps employed 
for the solution of equations (40). During the process four constants of inte- 
gration are introduced, namely, d-*\ <4'\ c','', c'/', and four are uniquely determined 
by the following conditions, the first of which determines c'^'~^' ; 

Ci.?o = , 2cy-» = ^'.*>„ - 2ZK.')„ , 

2c<;-» = ^,%-|-2(7^->„, 3c</)=-2^(%. 

The solution of equations (45) is 

Pi = «(i?o + 4*' cos T + (4*> sin T + • • • 

+ dg% cos {g + hm)r + ^;>^ sin (gr + Am)T + • • • , 

(47) 

w^ = c</^ + S[^l sin t + 7^*^ cos t -|- • • • 

+ ay^sin(5' + hm)T + y(i)^cos{g + hm)T + •••, 

where the coefficients are given by the formulas : 

"0,0 t ti,0^ 

'^^'"-l-ig + hmfy^^'" g + hmj' 

'^^-^ -1 - (gr + Am)^\,^''* +^ + AwiJ' 

(48) 8\-)„ = - 2c(') - iAi:l + IZ^,;) , 

1 / 7><> \ 

"•" g + hm\g + hm^'^'^''')' 

This completes the proof of the statement that the solution may be con- 
structed to any desired degree of accuracy. Of the four constants of integra- 
tion which occur in the calculation of the coefficients of fi*, one, cf^'\ is deter- 
mined in the step in which it enters ; two more, c'^*' and d^\ are determined in 
the next following step; while the fourth, cf^\ is determined in the second 
following step. 

Case II. When a uniform integral exists. The method of constructing the 
solution in this case differs from the preceding in one particular. The general 
steps are the same and formulas (48) are applicable. The difference occurs in 
the determination of the constant c^^*~^^ In the proof of the existence of a 
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periodic solution of equations (19) it was shown that the origin of time is arbi- 
trary, which is equivalent to the statement that at t = t^ the value w = Wg may 
be assumed arbitrarily. It follows that (7[,'^o = ^ ^^^ <io6S iio* determine (//"'" . 
This constant is determined by assuming an initial condition, for example, w = 
at T = . Then c'*~^^ is determined in the step in which it enters by the con- 
dition Wj_j = at T = . When the step for computing the coefficients of /** is 
completed, there remain two constants of integration, (f^'> and (f*'' which have not 
been determined. These constants are determined in the next following step. 

Case I. Symmetrical solutions. In this case there are two arbitrary initial 
conditions and the method of constructing the solution differs from that of case 
III in two particulars. The differential equations are (11). Consider the first 
step in the construction [equations (37) and (38)]. We determine c^^ and c^'^ 
respectively by the conditions 

Therefore <fp = (fp = . And (fp is determined by the condition that w^ shall 
be periodic while c^j'' is not determined until the next step. 

This process is applicable to all the following steps. The differential equa- 
tions (11) have a particular form which admits a symmetrical solution, and it 
may be established by complete induction that the equations for the determina- 
tion of p. and M)j have the form 

^-2^1 -3p,=^<-)o-H(V-4'-''+^i:'o)coST+ . . . + A^j^,cos{g+hm)T+. ■ ■, 
(49) 

where c'^* ~'' is determined by the condition 

The solution of (49) is 

Pi = o'oM) + <^i'' cos T 4- \- a!^%coa (g + hm)T + . ■ ■ , 

■Wi = ^i'.i) sin T -I- h ^j,<\ sin {g + hm)T + ■■ ■, 

Mi) _ _ JL /!(•) 

"o, — J-'^O, 5 

1 / 2i>*\ \ 

'^^.*-l_(^ + Am)^V *■" g+hmj' 

Si;)=-24^'-f^(;)-hJi)i;i, 

*"■*- g + hm\g + hm^'^'^'-'')- 



where 
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Numerical Examples. 

As illustrations of the preceding analysis three numerical examples will be 
treated briefly. It has not been shown that the processes are valid for the 
numerical values which are employed and which have been selected for conven- 
ience in graphical representation. It is probable that the series are convergent, 
although it has not been found possible to determine the true radii of convergence. 

In each case the mass of Jf is taken as the unit of mass, and M^, of mass 10, 
is supposed to move about M ma, circle whose radius is taken as the unit of 
distance. The unit of time is chosen so that the angular velocity N of Jf, is 
unity. The orbits considered are those in which the particle revolves about M. 
The period is chosen so that v = !y, whence fi = 0.2 . 

In the first example M^ of mass 5 is supposed to revolve about Jf in a circle 
so that the three finite masses remain always in a straight line in the order 
M^MM^. It follows from the solution of the quintic equation of Lagkange 
that the distance M^M\% 0.77172 • • • . The orbit of the particle is symmetrical 
with respect to the line joining the bodies. 

In the second example M^ of mass 6 is supposed to revolve about Jf in a 
circle so that the three bodies remain always at the vertices of an equilateral 
triangle. The differential equations admit the integral of Jacobi, but there is 
no symmetry theorem. 

In the third case it is supposed that M^ of unit mass revolves about J/j in a 
circle of radius A^. The orbit of the particle is symmetrical with respect to the 
radius vector of Jf, and crosses this line when the finite bodies are in conjunction. 

The differential equations of motion. With reference to Jf as origin and 
an axis passing always through Jfj let the coordinates of Jf, , M^ , and P be 
respectively ( 1 , ) , ( -Bj , W^ ) , and {r, w). The equations of motion, see equa- 
tions (7), of the particle are 

<pp ( dw X" 1 

lOi — P\ir + f^) + -2 = ^''J^,/iV[J{l+3cos2t«}4-|ap{3cosw-f6cos3«>} 

+ T6«V { 9 -H 20 cos 2w -I- 35 cos 4w } -f . . .] 
+-j^ mV[J{1+3 cos 2(«,- W,)} + l^p{Z cos {w- W^)+b cos Z{w - W,)} 

+ T^5 f;^) pM 9 + 20cos 2(m> - TF,) -f 35 cos 4(r« - W,)] -I- • ■ 1 , 

p --J-; + 2t-(t — |-m)=— I^M^ m^p [ J sin 2w + ^ap { sin m) -f- 5 sin 3w> } 

+ -^^a^p^ { 2 sin 2m) -H 7 sin 4w } -|- • • •] 

/fc^jf r a 

~ ~lif '''''' [^ sin 2(m,- TTJ + f ^P (sin {w - TTJ + 5 sin3(«r- TF,)} 
+ I'^ff ( J-JVm 2 sin 2(«, - TFJ 4- 7 sin 4( «, - TTJ} -f . . . j. 
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where a is given by tlie relation i^^ = Ji?. The numerical value of Jt? depends 
upon the masses and orbits of the finite bodies and must be determined sepa- 
rately for each of the three examples. 

Example 1. In this case R^ = 0.77172 and W.^ = tt. Now 1^ is given by 
the relation 



^,_ M,-\-M^M^ 



(^+(-r^)^^ 



Jf + M, + M,R, 
whence 

A;2 = 0.23763. 
It follows that 

PJfj = 2.37630 , -™-* = 2.58618 , 



R 



2 



o = 1.06914/i, ^ = 1.37222/*. 

-«2 






The equations of motion become 

' -pi^j^S +\ = (2.48074 + 7.44222 cos 'hw)py? 

— (1.15938 cos w + 1.93230 cos 3w)/3V* 

+ (4.23765 + 9.41700 cos 2w + 16.47975 cos 4«> )/>>* + •• -, 

*'°^S? + 4,(*+ ") — ('-"^^^ - ^'°>''' 

+ (0.38646 sin w + 1.93230 sin 3m»)/j>« 

- (4.70850 sin 2w +. 16.47975 sin 4m>)/j»/** + • • • • 

The periodic solution of equations (50) is 

/3 = 1 - I/* - (0.27136 + 0.96615 cos t + 4.96148 cos 2t)(i' 

+(0.62584-19.13740 cost-2.47963 cos 2t+ 0.40256 cos3tW+ • • •, 
^ ' w = T + (1.93230 sinT + 6.82204 sin 2t)m' 

+ (41.10885 sin t + 10.74793 sin 2t - 0.48307 sin 3t>» + • ■ • . 

Example 2. In this case -Bjj = 1 and W^ = irjZ . Here ¥ is given by the 
relation 

whence 

F = 0.06250. 

It follows that 

Fi!f, = 0.62500, PJ!fj = 1.56260/*, a = 0.67860/*. 



dr'' 
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The equations of motion become 

'-/,(^ + ^y+i-= (0.31250 + 0.93750 cos 2w)p/** 

4- (0.78125 - 1.17188 cos 2w + 2.02977 sin 2w)/3/*» 

+ (0.47714 cos w + 0.79523 cos 3to)/t>V* 

+ (0.59643 cos «> + 1.03308 sin w — 1.98808 cos 3w)pV* 

+ (0.16188 + 0.35974 cos 2w + 0.62954 cos 4w )/»>*+•• •, 
(62) ^ 

/>^ + 2|(^ + M) = -(0.937508in2«,)p^^ 

+ (1.17188 sin 2«> + 2.02977 cos 2w)/>/*» 

- (0.15905 sin w + 0.79523 sin 3w)/3V* 

+ (— 0.19881 sin w + 0.34436 cos w + 1.98808 sin 3«>)p>* 

- (0.17987 sin 2to + 0.62954 sin 4w)/3V' + • • • 

The periodic solution of equations (52) is 

p = 1 - f^ + (0.45139 + 0.89762 cos t - 0.62500 cos 2t)m' 

+ ( —0.47647+1.04542 cos t + 0.86090 sin t+0.46875 cos 2t 

- 1.36318 sin 2t - 0.16567 cos 3t) /*» + •• •, 
^ ^ tt, = T + (- 0.79524 sin t + 0.85938 sin 2t)/*« 

+ (0.13882 - 3.25720 sin t + 1.72180 cos t + 0.27995 sin 2t 

- 1.86062 cos 2t + 0.19881 sin 3t )/*» + •• •. 

Example 3. In this case M^ revolves about Jf, in a circle of radius A^ with 
an angular velocity N^. Let N^ be determined by the relation 

A^, jBj,, and W^ are determined by the relations 



B^ = 1/1+ A\ + 2^jCos2v< =\/1+A\+ 2^, cos 2t, 

A, sin 2t 



,;„Tr,= -'-^ 



B\ + l-A\ 



cosF;= 2iB, 



186 W. R. LONGLET : PERIODIC ORBITS OF AN INPINITE8IMAL BODY [April 

Here i? is given by the relation 

whence 

F = 0.09091 . 
It follows that 

PJfi = 0.90909, Pjtf, = 0.09091, 

A^ = 1.01200/t, a = 0.76910yit. 

The differential equations become 






(dw \* 1 
-3Z + f^) + -2 = (0.50000 + 1.50000 cos 2w)pn^ 



+ (0.86523 cos w + 1.44205 cos 3w)pV' 

+ (0.33273 + 0.73940 cos %w + 1.29395 cos 4w)/jV* 

+ (0.27600 sin 2t sin 2w — 0.13800 cos 2t — 0.41400 cos 2t cos 2«? )/»/*» 

+ (0.10474 + 0.17456 cos 4t + 0.31422 cos 1w + 0.10471 cos 4t cos 2to 

— 0.65864 sin 4t sin 2w)piJi* 

+ (0.07960 sin 2t sin w — 0.31840 cos 2t cos w — 0.53068 cos 2t cos Sw 

(54) + 0.39801 sin 2t sin 3w)pV* + • • • , 

d^w ^dpfdw \ ,^ ,„„«« . « s 

^ ^ + ^ ^ (57 + '*) = -( ^-^^^^^ '''^ 2"')'''* 

- (0.28841 sin w + 1.44205 sin 3«))pV' 

- (0.36970 sin 2w + 1.29395 sin 4w)/>V* 
+ (0.27600 sin 2t cos 2m) + 0.41400 cos 2t sin 2w)/o/*' 

— ( 0.31422 sin 2w + 0.10471 cos 4t sin 1w + 0.56864 sin 4t cos 2w)pA'* 
+ (0.02663 sin 2t cos w + 0.10613 cos 2t sin w + 0.39801 sin 2t cos Zw 

+ 0.63068 cos 2t s,mZw)p^ p.* + ■■ .. 
The periodic solution of equations (54) is 

p = 1 - I/* + ( 0.38889 + 0.72102 cos t - cos 'ir)p? 

+ ( - 0.02616 + 2.09168 cos t - 0.45400 cos 2t - 0.30043 cos 3t 

+ 0.03450 cos 4t) /*» + ..•, 
(65) yf* -r . 

w = T + ( - 1.44204 sin t + 1.37500 sin 2t) /*=' 

+ ( - 6.29838 sin t + 2.12066 sin 2t + 0.36051 sin 3t 

-0.03881 sin 4T)/t» + .... 
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Fig. 2. 
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The orbits represented by the solutions (51), (53), and (56) are illustrated by 
figures 2, 3, and 4 respectively. The comparison circles in the figures are not the 
circular orbits which have been called the undisturbed orbits. The undisturbed 
orbits are referred to fixed axes while the drawings are made with reference to 




rotating axes. The comparison circles represent orbits in which the particle 
would make a complete revolution with respect to the rotating axes during the 
period. In figure 2 the points which are numbered 1, 2, • ■ • represent the posi- 
tions of the particle in the periodic orbit at intervals in t of 7r/12. The cor- 
responding positions in the comparison circle are indicated by the numbers 
1', 2', . • -. In figures 3 and 4 the intervals are 7r/4. 



